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AN EXPERIMENTAL INVESTIGATION IN AN ATMOSPHERE 
ENTRY SIMULATOR OF NYLON AS AN ABLATIVE 
MATERIAL FOR BALLISTIC MISSILES* 

By Gary H. Bowman and Raymond C. Savin 


SUMMARY 


Tests were conducted in a small-scale atmosphere- entry simulator to 
investigate the feasibility of nylon as an ablative material for long- 
range ballistic missiles. The test models were spherically tipped cone 
cylinders 0.22 inch in diameter. The cone half-angles varied from 40° 
to 90 the tip radius was equal to approximately one fourth of the 

model diameter. The tests simulated the motion and heating associated 
with the atmosphere entry of ballistic missiles having entrance veloci- 
ties ranging from 13,300 to 19,800 feet per second, diameters from 2.6 
to 2.9 feet, weights from 707 to 930 pounds, and maximum, ranges from 
12l8 to 3^00 statute miles, respectively. 

Studies of models recovered after flight through the simulator 
indicate that nylon is capable of withstanding the thermal stresses 
associated with atmosphere-entry heating. These studies also indicate 
that increasing cone angle causes, for the most part, a slight decrease 
in ablation; whereas, increasing entry velocity causes an increase in 
ablation. The weight loss due to ablation is about 2 percent or less of 
the initial model weight in all cases tested, which indicates that nylon 
is an effective ablative material for long-range ballistic missiles. 


INTRODUCTION 


The aerodynamic heating associated with the atmosphere entry of a 
ballistic missile can impose severe weight penalties on a vehicle since 
some method of structural cooling must be employed. Currently one of 
the more promising methods of providing the required cooling is the use 
of an ablation heat shield. With this method the forward face of a 
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missile is covered with expendable material which, as it ablates, absorbs 
heat by essentially utilizing its heat capacity and latent heat of vapor- 
ization. In addition, the vapor injection into the air boundary layer 
may also block a large fraction of heat from reaching the surface. 

Recent theoretical and experimental investigations (see, e.g. , 
refs. 1 and 2) have indicated that some materials belonging to the 
plastics group have attractive ablation characteristics. Unfortunately, 
the ablation cooling provided by many of these materials (e.g., nylon) 
is as yet not amenable to accurate theoretical treatment since applica- 
tion of theory depends to such a great extent on a knowledge of unavail- 
able material properties in various states (i.e., solid, liquid, and 
gaseous). Although there are some experimental data available for nylon 
and other materials, these data have been obtained, in general, under 
conditions which do not simulate the history of the heating associated 
with the atmosphere entry of a full-scale vehicle. As a result, certain 
important aspects of material capability (e.g., structural ability to 
withstand thermal stress due to this heating) have in most cases not been 
determined. Another example would be the net weight loss over the 
trajectory. 

The motion and heating of a ballistic missile during atmosphere 
entry can be simulated with the small-scale atmosphere-entry simulator 
at the Ames Research Center. Thus, the simulator is particularly suited 
for studies in which heating history plays an important role. Accord- 
ingly, an exploratory investigation of the effectiveness of nylon as an 
ablative material for long-range ballistic missiles was conducted with 
this equipment. The results of this investigation are the subject of 
this report. 


NOTATION 


Cp specific heat at constant pressure 

Cg specific heat of material up to melting point 

Cp drag coefficient referred to base area 

D diameter of model or missile 

h specific enthalpy 

h l heat capacity of solid-liquid system 
k thermal conductivity 

l length 
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latent heat of fusion 

latent heat of vaporization 

mass rate of ablation per unit area 

mass rate of vaporization per unit area 

molecular weight ratio of air to material vapor 

static pressure 

Prandtl number 

local heating rate per unit area 
average heating rate per unit area 
total convective heat absorbed 
effective heat of ablation per unit mass 

average effective heat of ablation per unit mass over trajectory 
nose radius of model or missile 
missile range 

free-stream Reynolds number referred to hemispherical, nose radius 
of model 

distance along gas-liquid interface 
nose surface area of model or missile 
time 

temperature 
melting temperature 
velocity 

average rate of ablation weight loss per unit area 

initial weight of model or missile 

weight lost due to ablation 

distance along axis of model or missile 
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distance normal to gas-liquid interface 
thermal diffusivity 

conical half-angle of model or missile 

inclination of flight path from horizontal 

viscosity 

density 

shear stress 

ratio of heat transfer with and without vaporization 

Subscripts 

entrance conditions for simulated and real atmosphere 

conditions at gas-liquid interface 

material property in liquid state 

property associated with missile 

property associated with model 

nonablating surface 

material property in solid state 

stagnation conditions 

material property in vapor state 

APPARATUS AND TESTS 
Test Apparatus 


Tests were conducted in the Ames small-scale atmosphere-entry 
simulator. A detailed description of this equipment may be found in 
reference 3- In the present tests, the single-stage light-gas gun was 
employed to obtain model velocities (relative to the air stream) up to 
approximately 16,000 feet per second. To obtain higher velocities 
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(up to 20,000 feet per second), the gun was modified to in clude a 
second-stage pump tube and, thereby, provide a means of compressing the 
driving gas, helium, to higher pressures and temperatures. 


Models 


The test models were 0.22 caliber and were machined from nylon 
(Polypenco 101) . They consisted of a spherically tipped conical nose 
and a cylindrical afterbody. The cone half-angles employed were 4o°, 
50°, 60°, 70°, and 90° • The nose radius was held constant and equal to 
0.053 inch or approximately one fourth of the afterbody diameter. Since 
the forebody dimensions changed with cone angle, the afterbody length 
was varied to maintain the same weight of approximately 0.105 gram for 
all models. Dimensions of the models are given in figure 1. 


Test Procedure 


Prior to each test, the model weight and afterbody length and 
diameter were measured. All models were launched and a time-distance 
history was determined in the manner described in reference 3 • Entrance 
velocities ranged from approximately 13,300 to 19,800 feet per second 
relative to the air stream. The terminal velocity was maintained nearly 
constant at approximately 1000 feet per second. The models were recovered 
at the end of their flight in soft balsawood catcher material. The 
recovered models were then weighed and measured to determine their weight 
losses and any changes in dimensions. 


Reduction of Data 


There are three possible causes of weight loss from a model fired 
in simulation: (l) launching, (2) recovery, and (3) flight through the 

simulated atmosphere. The latter or ablation weight loss can be deter- 
mined only as a result of a careful evaluation of the weight losses due 
to launching and recovery and subtraction of these losses from the total 
measured weight loss. The weight loss due to launching was determined 
from the changes in afterbody length and the decrement in diameter due 
to the scraping of the model along the launch tube. The change in length 
caused by the hot helium gases scrubbing the base of the model during 
launching was found to be negligible. The weight loss due to recovery 
of the model was determined by a series of penetration tests. Models of 
the same weight and shape were launched at velocities comparable to the 
impact velocities of the simulation tests (approximately 1000 feet per 
second) and allowed to penetrate soft balsawood blocks - the same type 
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of material employed as catcher blocks for the tests. These tests 
indicated a negligible weight loss due to recovery. Accordingly, the 
weight lost by ablation during atmosphere entry was determined by sub- 
tracting the weight loss due to launching from the total weight loss of 
the test model. In general, the ablation weight loss amounted to about 
60 percent of the total weight loss. 


Accuracy of Test Results 


Errors in determining model weight losses can arise principally 
from the following two sources: (l) the measurement of the afterbody 

length and diameter and (2) the weighing of the models. The afterbody 
length of each model was measured to ±0.0002 inch prior to its flight 
through the simulator and the diameter was measured to ±0.0002 inch 
before and after the test. Each model was also weighed to ±0.00005 gram 
before and after the test. These limits of measurement result in a 
maximum probable error of approximately ±20 percent of the average 
minimum ablation weight loss observed. 


RESULTS AND DISCUSSION 


It has been demonstrated in reference 4 that the total convective 
heat per -unit mass and resulting thermal stresses associated with the 
atmosphere entry of a ballistic missile can be duplicated with a small- 
scale test model in an atmosphere-entry simulator. The simulation 
requires that the model and missile surfaces be geometrically similar 
and made of the same material and that the velocities and local Reynolds 
numbers be the same at corresponding points in the trajectory. The degree 
to which the present test conditions simulate the motion of full-scale 
ballistic missiles is indicated in figure 2. In this figure, theoreti- 
cal velocities for a simulated missile along with experimental velocities 
obtained for the test model are shown as a function of simulated altitude. 
The simulation of altitude requires the product of ambient density and 
vehicle diameter (i.e., pD) to be the same for model and missile. The 
results shown in this figure are representative of the test results in 
general. The data points were determined from one time-distance history 
obtained for the 50° cone half-angle model launched at a velocity (rela- 
tive to the air stream) of 15,600 feet per second. The dimensions of 
the simulated missile shown in the figure were determined according to 
reference 4. The entry angle, 9 %, and corresponding range were determined 
from the method of reference 5 and the theoretical velocities were calcu- 
lated according to reference 6. The agreement between the theoretical 
and experimental velocities indicates that the requirement of equal 
velocity for model and missile at corresponding points in the trajectories 
is essentially satisfied. The experimental velocities are somewhat lower 
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than the theoretical values near the end of the trajectory. However, 
the terminal velocities are small in both cases. Since the total heat 
absorbed per unit mass is proportional to the difference in the square 
of entrance and terminal velocities (see ref. 4), the error in simulat- 
ing the total heating per unit mass is negligibly small for our purposes. 

The range of the present investigation is indicated in the following 
table which lists the characteristics of the missiles whose atmosphere 
entry was simulated. 


&C, 

deg 

v E . 

ft/sec 

0 E> 

deg 

Diameter, 

ft 

Weight, 

lb 

Range, 

statute miles 

40 

15, 

600 

38 

2 

•7 

760 

I806 

50 

60 

70 

90 

\ 

t 

\ 

f 

\ 

> 

\ 

f 

\ 

' 

50 

13, 

300 

4i 

2, 

.6 

707 

1218 


17, 

300 

37 

2.8 

805 

2318 


19,800 

33 

2.9 

930 

3400 


These characteristics were determined in the manner just described. 

Thus, for the first group listed, the tests simulated the flight of 
missiles having shapes similar to the test models with a diameter of 
2.7 feet, a weight of 760 pounds and, according to reference 5, a total 
range of 1806 statute miles when entering the earth's atmosphere at 38° 
to the horizontal at a velocity of 15,600 feet per second. For an 
entrance velocity of 19,800 feet per second, the missile will be slightly 
larger, will weigh 930 pounds, and will have a range of 3400 statute 
miles . 

Let us first consider the results of the tests from the standpoint 
of the structural feasibility of nylon as an ablation heat shield for 
ballistic missiles. One qualitative measure of this feasibility is 
readily obtained from an examination of the test models after their 
flight through the simulator. In this connection, the majority of models 
recovered showed no evidence of cracking even at test velocities of 
19,800 feet per second. Because of the infrequency of cracking, it is 
believed that the cracks were caused by the stresses imposed in launching 
the models. 1 Aside from the occasional cracking, the models also experi- 
enced surface pitting. As indicated by the photographs of the 50° cone 
half-angle test models in figure 3, the pitting was more severe at the 
higher velocities. Conclusive evidence as to the cause of this pitting 

■^he test models experienced accelerations up to approximately 
15 million g's during launching. The decelerations experienced by the 
models, being of the same order as the accelerations, were also higher 
than those associated with a full-scale missile. 
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is lacking at present. More study is needed to determine whether it was 
caused by surface spalling or, as noted in reference 3 with copper-faced 
models, by impurities in the air stream. 

It should be mentioned that the solid models employed in the tests, 
being more rigid, will tend to fail more readily under thermal loading 
than would a relatively thin shield of the same material on the face of 
the model. In this connection, the analysis presented in appendix A 
indicates that the depth to which heat was conducted into the models was 
much less than 1 percent of the model diameter. This depth (a fraction 
of an inch for the simulated full-scale missile) is considerably less 
than might be required for a practical heat shield. Thus, the presence 
of the additional ablation material in the core of the model did not 
affect the heat capacity of the models in the tests. It would appear, 
then, on the basis of these qualitative results, that nylon is an effec- 
tive heat shield material. 

The next point to be considered is the amount of nylon required as 
an ablation heat-sink material. An indication of the amount of surface 
material lost due to ablation is shown in figure k. These photographs 
were made by projecting the silhouette of recovered models (S c = 50°) 
onto a screen containing a template of the pretest model profile . Since 
the proper orientation of the template with respect to the recovered 
models was difficult to achieve, the photographs do not provide an accu- 
rate measure of the total amount of material lost. They do indicate, 
however, that the greatest thickness of material was lost from the spheri- 
cal segment of the nose . It appears also that ablation has taken place 
over a greater area of the nose in the case of the model launched at the 
higher velocity. It is further indicated in this latter case that nose 
material has been lost in a localized region near the sphere-cone juncture. 
An ablation distribution of this type might be expected in the presence 
of transition to turbulent boundary-layer flow, since for this type of 
flow the heat-transfer rate reaches its maximum value downstream of the 
stagnation point (see ref. 7)* It was not possible to determine the type 
of boundary layer on the models during flight. However, it may be deduced 
from the free- stream Reynolds number values shown in figure 5 that turbu- 
lent flow might well have existed in the region of the sphere-cone 
juncture at the higher entrance velocities. 

The actual measured weight losses due to ablation of the test models 
are shown in figures 6 and J. The variation of percent weight loss with 
model shape is shown in figure 6. The drag coefficients indicated on the 
figure were computed from Newtonian theory. Each datum point represents 
the average weight loss obtained in launching four to five models at an 
entrance velocity of about 15,600 feet per second. The root-mean- square 
deviation from the average is also indicated. There are two things to be 
noted from this figure. First, one may speculate that the minimum weight 
losses observed indicate that the 60° and 70° models approach shapes having 
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minimum heat transfer. Perhaps the point of greater importance, however, 
is that the ablation weight loss appears to be less than 2 percent of the 
initial model weight for a rather large range of nose shapes. The impli- 
cation of this result is that one has a rather wide latitude in choice of 
shape at a cost of a relatively low weight penalty. The effect of 
entrance velocity on the ablation weight loss is shown in figure 7* These 
data were obtained as a result of launching 50 ° cone half-angle models 
over a range of velocities from 13,300 to 19,800 feet per second as indi- 
cated. Again the weight loss is about 2 percent or less over the velocity 
range. On the basis of these and previously discussed results, nylon 
appears to be an effective ablation heat shield material. 

It is appropriate now to consider the test results in terms of the 
effective heat capacity of nylon. To this end, consider an effective 
heat of ablation parameter, Q*, defined as the ratio of heating rate per 
unit area at the nose of a nonablating body to the rate of weight loss 
per unit area. Then, in terms of time-averaged quantities. 



w 


where the average heating rate is simply 



and Q is the total convective heat transferred to the nose area, S, 
over the trajectory time, t, and 


w 


AW 

ts 


where AW is the total weight lost due to ablation. Therefore, Q* can 
also be expressed in the form 


AW/W 

Accordingly, if Q* mo = Q*mi> the fraction of weight lost will be the 
same in model and missile since the total heat per unit weight is 
duplicated (ref. 4). 

The data for the 50° test models have been reduced to values of Q* 
according to the above relation. The values of Q, the total heat trans- 
ferred to the model, were calculated for two cases: (l) the flow 
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downstream of the stagnation point was all laminar and (2) the flow was 
all turbulent. The laminar heating rates were calculated in the manner 
of reference 8 and the turbulent rates as suggested in reference 7- The 
stagnation heat flux was determined according to the analysis of refer- 
ence 9 . The transient time required for the onset of ablation was neg- 
lected since the analysis of appendix A demonstrates this approximation 
to be justified. The variation of average effective heats of ablation 
for both laminar and turbulent heating rates is shown plotted in figure 8 
as a function of the stagnation enthalpy, h-fcgj corresponding to the 
initial point in the trajectory. Because of the relatively high veloci- 
ties under consideration, the relation 



was employed in determining the specific enthalpies. Since the entrance 
velocity is the same for model and missile, the stagnation enthalpy is 
essentially duplicated. 

An analysis of the degree to which simulation of the ablation process 
can be achieved in tests of the present type is presented in appendix B. 

It is indicated by this analysis that if the ablation process is purely 
one of liquefaction then the test results for Q* are too large and must 
be corrected before they can be applied to a full-scale missile; however, 
if ablation takes place solely as a result of vaporization, then the test 
results are essentially applicable without correction. In order to esti- 
mate the correction in the case of liquefaction, it is necessary to know 
the thermal properties of the ablating material in its solid and liquid 
states. Unfortunately, there appears to be a dearth of information con- 
cerning the liquid-state properties, not only for nylon, but for thermo- 
plastics in general. Calculations were carried out for pyrex glass for 
which the pertinent properties are available. The results are shown in 
figure 9, where a maximum correction of about 10 percent in Q* is indi- 
cated for scale factors corresponding to the present tests (i.e., 

D mj /Dmn ~ 150 ) . It should be noted that nylon has a melting temperature 
of about 500° to 600° F. Some data presented in reference 10 indicate 
that the vaporization temperature of nylon, or perhaps more properly 
depolymerization temperature, is approximately 600° to 700° F in a vacuum. 
Although one would expect higher depolymerization temperatures at higher 
pressures, it is believed that the surface could not maintain much higher 
temperatures without depolymerizing. In this connection, it should be 
mentioned that little, if any, evidence of melting could be detected from 
recovered test models. It is reasoned on this basis that the depolymeri- 
zation temperature is close to the melting temperature and therefore 
essentially all the ablation will take place by surface depolymerization 
in both model and missile. Thus, for laminar flow, the conclusion of 
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case (b) in appendix B 2 (i.e., my/riiL -*• l) is applicable to the present 
test results; 3 namely, Q* mo ~ Q* mi . In the case of turbulent flow, the 
situation may not be the same. However, it does not appear illogical to 
assume comparable simulation (to the laminar flow case) in view of the 
fact that the local Reynolds numbers and driving enthalpies are essen- 
tially duplicated in model and missile. 


CONCLUDING REMARKS 


Tests have been conducted in the small-scale atmosphere- entry 
simulator to assess the feasibility of nylon as an ablation heat shield 
material. The test models consisted of spherically tipped cone cylinders 
which had a diameter of 0.22 inch and weighed approximately 0.1 gram. 

The cone half-angle varied from ^0° to 90° aiii the tip radius was equal 
to approximately one fourth the afterbody diameter. The tests simulated 
the motion and heating associated with the atmosphere entry of ballistic 
missiles having entrance velocities ranging from 13,300 to 19,800 feet 
per second, diameters from 2.6 to 2.9 feet, weights from 707 to 930 pounds, 
and maximum ranges from 1218 to 3^00 statute miles, respectively. 

Nylon was found to be a structurally attractive material from the 
standpoint of withstanding the thermal loading associated with the atmos- 
phere entry of ballistic missiles. Ablation weight loss was found to 
decrease with increasing cone half-angle up to 60° and remained essen- 
tially constant or increased slightly for blunter shapes. Increasing 
velocity caused an increase in ablation. The measured weight loss due 
to ablation was about 2 percent or less of the initial model weight in 
all cases. On the basis of a qualitative simulation analysis also pre- 
sented, it was indicated that these results are approximately applicable 
to full-scale ballistic missiles. In view of these findings, it is con- 
cluded that nylon is an effective ablation heat shield material for 
long-range ballistic missiles. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., Aug. 13, 1959 


2 The conclusions of appendix B for the stagnation point can be shown 
to hold downstream along the body surface by means of Lees' "local simi- 
larity” assumptions (ref. 8) which have been, in essence, verified 
experimentally (ref. 11). 

3 It should be noted from appendix A that the transient time is 
simulated when Q* mo = Q* m i- 
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APPENDIX A 

ESTIMATE OF THERMAL THICKNESS AND TRANSIENT 
TIME FOR ONSET OF ABLATION 
Initial Surface Temperature Rise 


The variation of surface temperature with time, t, on a semi- inf inite 
solid may he expressed in the form (see ref. 12) 


T 




(Al) 


where k is the thermal conductivity and a is the thermal diffusivity 
of the solid. Now the heat-transfer rate, q Q , is larger in the simulator 
hy the ratio of missile to model diameter (see ref. 4). If we require the 
surface temperature to be the same in model and missile, then according 
to equation (Al) , t must be smaller in the simulator by the square of 
the ratio of model to missile diameter. This relation is identical to 
the simulation requirement for the trajectory times of model and missile 
obtained in reference 4. Therefore, prior to the onset of ablation, the 
fraction of trajectory time required for the surface to heat up is 
duplicated in model and missile. 

The variation of surface temperature at the stagnation point with 
time for the trajectory corresponding to Vj; = 15,600 feet per second 
was calculated from equation (Al) . The following thermal properties of 
nylon were employed in these and all subsequent calculations: 
k = 0.1405 Btu/hr-ft °F, c p = 0.4o Btu/lb °F, p = 68.6 lb/cu ft. Since 
the heat-transfer rate, q Q , varies with time, the integral was evaluated 
graphically. The results are presented in figure 10 where it is evident 
that the surface temperatures become large in a matter of microseconds. 

If the average heat-transfer rate over the nose had been used to calcu- 
late the surface temperatures, the resulting times would be approximately 
four times those shown. They would still be negligible, however, compared 
to the test flight *time of 1 millisecond. 


Thickness of Thermal Layer 


Since the initial temperature rise occurs over such a small interval 
of time, the heating rate can be assumed constant during this period. 

The relation defining the temperature at each point in a homogeneous 
semi- infinite solid at time t is then given by (see ref. 12) 
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T 



4at 


- x erfc 


x 

2^/at 


(A2) 


It should be noted that relations (Al) and (A2) become identical for 
x = 0 and q Q = constant. The temperature distribution through the model 
was calculated from equation (A2) for t corresponding to a surface 
temperature of 1000° F (see fig. 10) . The stagnation heating rate cor- 
responding to the initial point in the trajectory was used in the calcu- 
lations; namely , at V = V E = 15,600 ft/sec, (4 0 ) t = 35>000 Btu/sq ft-sec. 
The results are shown in figure 11. For convenience, the initial temper- 
ature of the model interior has been taken as 0° F. It is- clear from 
these results that the heat conduction effects are confined to a region 
of thickness approximately 0.02 percent of the model diameter. This 
percentage thickness is substantially less than the shell thickness 
required to simulate a practical heat shield for a full-scale missile. 


Transient Time for Steady-State Ablation 


After the surface has reached a temperature sufficiently high for 
the onset of ablation, there is a further time lag or transient time for 
ablation to attain a state of equilibrium. An order of magnitude esti- 
mate of this transient time may be obtained from the relation (see 
ref. 13) 


t 


k s p s 

c s 



(A3) 


It is interesting to note from this expression that if Q* mo = Q* m i> the 
transient time will be smaller for the model by the square of the ratio 
of model to missile diameters (note q Q ~ l/D from ref. 4). Therefore, 
the fraction of trajectory time required for the transient time will be 
the same for model and missile. A typical transient time for the present 
tests has been calculated employing an average value of q 0 and for the 
trajectory corresponding to &c = 90° and Vg = 15,600 ft/sec; namely, 

4 0 = 12,000 Btu/sq ft-sec and Q* = 1000 Btu/lb (see fig. 8). The result- 
ing transient time is about 0.05 millisecond which is about 5 percent of 
the total flight time. 
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APPENDIX B 


THEORETICAL ESTIMATE OF ABLATION 
WEIGHT LOSS SIMULATION 


Relatively simplified expressions for steady-state surface ablation 
phenomena have recently been obtained by Lees (ref. 7) and. Bethe and 
Adams (ref. l4) . The purpose of the following analysis is to utilize 
the results of these studies in an effort to obtain an estimate of the 
degree to which simulator test results are applicable to full-scale 
missiles. For details of the development of the basic equations and for 
the simplifying assumptions involved, the reader is referred to the refer- 
ences cited above. For simplicity, the following analysis is restricted 
to the stagnation region of a blunt body. 

With these points in mind, and on the basis of the analyses of 
references 7 and l4, the energy, momentum, and continuity equations may 
be combined to yield 


m L - my 


2 2 * 2 . 
A C P L ^Li m L 


(Bl) 


This expression states that the mass rate of ablation less the mass rate 
of vaporization (i.e., the quantity of liquid flowing off the nose) is 
proportional to the liquid density, the square of the average thermal 
conductivity of the liquid, the square of the temperature at the gas- 
liquid interface, the shear gradient at the interface, and is inversely 
proportional to the liquid viscosity at the interface, the square of the 
heat capacity of the liquid and the square of the mass-rate of liquefac- 
tion. The coordinate system to which this and subsequent equations refer 
is stationary with respect to the gas-liquid interface as shown in the 
following sketch. 
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In the presence of vaporization, the aerodynamic heat-transfer rate 
at the interface must be balanced by the heat of vaporization in addition 
to the heat capacity of the solid and the liquid. Therefore, the energy 
balance equation may be written 

4i = ^l h L + 

where the heat capacity of the solid and liquid is 

H L = L S + c s (T ra - T S ) + Cp L 

and Lg is the latent heat of fusion and Ly is the latent heat of 
vaporization. 

The effect of mass injection into the boundary layer is to reduce 
the heat-transfer rate and shear at the interface. Therefore, we can 
write 


(B2) 

(B3) 


ii 

io 



(b4) 
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where \|r is often referred to as the "blowing" parameter and has been 
the subject of numerous investigations (see, e.g., refs. 15 through l8) . 
This parameter may be defined in terms of the driving enthalpy and mass 
injection rate as (see ref. l4) 


t = l-aM b (h t -h i ) (B5) 

To 

where M is the ratio of the molecular weights of air to the injected 
vapor and a and b are empirically determined constants. The rate of 
mass injection may be related to the vapor pressure by (ref. l4) 


^ = l-aM*" 1 + 
PV 


To 


M(h t -h i )m v 


(B6) 


At the stagnation point the shear gradient may be written (refs. l4 
and 19 ) 


T o 


ht-hj. r 



Substitution of this relation along with expression (B4) into equation (Bl) 
yields 

! _ my = Zp&L (TjfMvtM Pr t 2/3 tqo (B?) 

m L ( H - hi ) rp Li m L 3 Cp L 2 A 2 

The combination of equations (B2), (b 4), and (B5) results in 


^2 = h l + ~ [L V + aM b (h t - hi)] (B8) 

m L ra L 

for the energy balance equation. Therefore, consistent with the definition 
of the effective heat capacity employed in the text 


Q* 


^ = H l + ^ [Ly + aM b (h t - h ± )j 
m L m L 


(B9) 


For materials whose viscosity and vapor pressure are steep functions 
of temperature we may write 
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PL = 




P V = 





where A , B, C , and D are constants. Then 


PL ~ PV 


-A/C 


and (as suggested in ref. l4) 


P Li 


Ai 


-A/C 


(BIO) 


(Bll) 


where is the viscosity at the temperature at which the vapor pres- 

sure is equal to the stagnation pressure. In view of equations (B9) 
and (Bll), expression (B7) finally may be written 


1 Q*~ H l = 2p L k L 2 (T i ) 2 j2(pt/pt) 

Ly+ aM b (h t - hj_) r(h t - h i )p Lt q 0 2 c 

where, by virtue of equations (b 6) and (B9), 

p t -j,! Q* ( ^ + . 

P Vi M Q* - H l W - Bq 


Pr t 2/3 Q* 3 ilr /p t \- A//C 


2 A 2 
P L A 

V p Vi/ 

(B12) 

aM b H L ^ 
Q* / 


(B13) 


and, from expressions (B5) and (B9)> 


aM b (h t -hi)(Q*-H L ) 
Q*[ Ly + aM b ( h t - hq ) ] 


(Bl4) 


Consider now equation (B12) in the following three limiting cases as 
it applies to test model and full-scale missile. 


Case (a) Liquefaction with no vaporization (my = 0): 

Since my = 0 for this case it follows from equations (B5) and (B9) 
that t = 1 and Q* = Hp and equation (B12) reduces to 


2 PL k L 2 ( T i^^ g ( p t/Pt) Pr t 2/3 Q* 3 

r(h t -h i )p Li q 0 2 Cp L 2 A 2 
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Now the density of the liquid, the stagnation Prandtl number, and the 
ratio of stagnation pressure to stagnation density will be the same for 
model and missile. Furthermore, the thermal conductivity of the liquid 
can be considered approximately equal in both cases since it varies 
relatively slowly with temperature compared to the viscosity. If we 
note that (ref. 4) 


( 4p ) mo 

(4o) m i 


Pmi 

Ano 


and, for large flight velocities. 


(h+ - hp) 
v _ - 1 mo 

(ht - hi) m p 


then equation (B15) may be expressed in the form 


Q mo _ ^^^mo 

1 mo 

1/3 


^ mi ( E L) m i 

(nil). 


(Tp) 

mo 


%i 

ftiO' 


1/3 


It is evident from equations (B3) and (BIO) that the above relation can 
be evaluated only for materials for which the thermal properties in its 
liquid state are known. For this reason, the variation of Q* mo /Q*mi 
with the scale factor D mi AVip for two values of the interfacial temper- 
ature for the missile have been calculated for pyrex glass. The results 
are shown in figure 9- The thermal properties of pyrex employed in this 
calculation were obtained from references 7 and l4. Thus, a melting 
temperature of 3000° R was employed. In figure 9, the curve Tp/T m =1.5 
corresponds to a flight altitude of 60,000 feet at a velocity of 17,000 
feet per second according to the calculations of reference 7- Based on 
an extrapolation of these calculations, the curve Tp/T m = 2 corresponds 
to a flight velocity in excess of 22,000 feet per second. Although the 
results shown in this figure are strictly applicable to pyrex glass, they 
should nevertheless hold qualitatively for materials which have a high 
liquid viscosity. From the standpoint of duplicating ablation by liquefac- 
tion of the surface material on a missile by means of a test model, the 
advantage of using large-scale facilities is evident from figure The 
missiles simulated by the present tests are approximately 150 times the 
size of the test models. It is indicated from figure 9, therefore, that 
the largest probable factor for the ratio of model to missile effective 
heat capacities is approximately 1.1. 
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Case (b) Strong vaporization with high viscosity (my/iriL -*■ l): 

If there is strong vaporization on both model and missile, then 

H l « Q* 

and expressions (3313 ) and (Bl4) reduce to 

In „ i + — h — 

P V ± M(h t -hj_) 


Ly+ aM b (h t - hp) 

respectively. Since the driving enthalpy, hp-hp, is essentially dupli- 
cated in model and missile, the above parameters will be approximately 
the same in both cases. 4 If, in addition, the liquid viscosity is high, 
then the right-hand member of equation (3312) will be very small and, in 
the limit, this equation may be expressed in the form 

q* [L V + aM b (h t - hi)] 

Q*mi [ Ly + aJVI b (h t - hp) ] m± 

Again the driving enthalpy is the important parameter. The result for 
Q* indicated by the above expression is similar to that obtained theo- 
retically in reference 14 and verified experimentally (for teflon) in 
reference 20. In this limiting case, therefore, Q* mo -*• Q*mi' 

Case (c) Strong vaporization with moderate viscosity (iiiy/mp < l): 

In this case, the term on the right in equation (B12) cannot be 
neglected. If the fraction of liquid vaporized is the same in model and 
missile, then 


and the left-hand member of equation (3312) is the same for both model and 
missile since from equation (339) 


4 Since the assumption of Lewis number equal to unity is employed 
(see refs. 7 and lh) , it follows that the molecular weight ratio, M, is 
equal in both model and missile. 
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Q*- H L = my 
Ly + aM b (h^- - hj_) 

Equation (B12) then becomes (for strong vaporization) 


(B16) 


( l- L Irt ) 

■ Lrt me 


1 / 3|— — l 2 / 3 

NJ / 


which, it will be noted, is similar to the limiting expression obtained 
for no vaporization (case (a)). In the present case, however, 


(d L +) 

^ mo 
( ^ } mi 


rwj 

|_^ p t^mJ 


so that 


Q mo 
Q*mi 


0 _ ^i^rni 

1 i T i)mo_ 


3(l-A/C) 


Since (Ti) mo /(Tj_) ^ > 1, it is clear that there is some minimum value of 
A/C (< l) for which Q* mo = Q* m i- For values of A/C greater than this 
minimum value , Q* mo < Q*mi and the test model results will be conserva- 
tive with respect to full-scale results. The term "conservative" is used 
in the sense that a greater percentage of material is lost on the model 
than would be lost on the simulated missile. In the case of pyrex glass, 
A/C ~ 0.80. For other materials for which A/C « 1, the test model 
results for Q* in the present case will be greater by a factor always 
less than that which was observed for no vaporization (case (a)). 


> /fV 


the left-hand member of equation (B12) will always be less in the case of 
the model. Therefore, the ratio Q*mo/Q*mi will be less than in the 
previous case just considered; namely (my/m^) = (my/m^) . . 
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Finally, for the case 'where 

Wao W4 l 

equation (B1 6 ) expresses the fact that (noting that Hp « Q*) Q* mo < Q*mi' 
Therefore the effective heat capacity determined from the tests will again 
be conservative with respect to the full-scale values. 

On the basis of the foregoing analysis for the three limiting cases, 
m v /m L = 0, m v /m L < 1, and my/iiiL ->-1, it is indicated that simulator test 
results are least applicable when no vaporization occurs. Thus, the 
results presented in figure 9 may he considered representative of an 
approximate maximum correction to simulator test results for Q* in order 
to apply them to the full-scale missile. 
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Figure 2.- Comparison of theoretical variation of velocity with altitude for full-scale missile 
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Figure 3 .- Typical preflight and postflight models; & c = 50 C 





ro 

^1 




CONFIDENTIAL 


ro 

oo 




• • i 

• • « 
• • • • « 


• ••*« 


« c 

I • • • • 


• • • • • 
• • 


*•• •• 
• • • • 


Unfired model 

Figure 4.- 
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Figure 6.- Variation of ablation weight loss with nose shape; V E = 15,600 feet per second. 



Figure 7-~ Variation of ablation weight loss with entrance velocity; 
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Figure 10.- Initial surface-temperature rise at the stagnation point of a 
nylon model during flight in the simulator; Vg = 15,600 feet per second. 
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Figure 11.- Temperature distribution in nylon model at beg inn ing of flight 
through simulator; V E = 15,600 feet per second. 
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